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ABSTRACT 

Data from the Lunar Exploration Neutron Detector (LEND) Collimated Sensors for Epither- 
mal Neutrons (CSETN) are used in conjunction with a model based on results from the Lunar 
Prospector mission to quantify the extent of the background in the LEND CSETN. A simple 
likelihood analysis implies that at least 90% of the lunar component of the LEND CSETN flux 
results from high energy epithermal neutrons passing through the walls of the collimator. Thus, 
the effective full-width at half-maximum of the LEND CSETN field of view is comparable with 
that of the omni-directional Lunar Prospector Neutron Spectrometer. The resulting map of 
high energy epithermal neutrons offers the opportunity to probe the hydrogen abundance at low 
latitudes, and provide constraints on the distribution of lunar water. 



Subject headings: Moon — methods: data analysis 

1. Introduction 

The weathering of the lunar surface is so 
slow that it should provide extractable informa- 
tion about the history of volatile molecule trans- 
port within the inner Solar System, a subject 
of considerable relevance to water-based life on 
Earth. Samples of lunar regolith returned in the 
Apollo missions contained a sufficiently low abun- 
dance of water that these molecules were, until 



recently, dismissed as terrestrial contamination 
([Epstein fc Tavlorlll973l) and the Moon was con- 
sidered to be dry. The change of paradigm has 
been driven by evidence that has accumulated 
over the past 15 years from a numbe r of differ 



ent missions and types of experiment ([Lawrence 



2011). 



The first of this recent wave of techniques for 
divining for lunar water was the use of bistatic 
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circularly polarised radar measurements from the 
Clementine mission (jNozette et alJ ll996). A peak 
at opposition in the ratio of returned same-sense to 
opposite sense circularly polarised flux was posited 
to be due to multiple scattering within a low-loss 
medium such as water ice. The signal correlated 
with Clementine being over the lunar poles, the 
locations of permanently shaded 'cold traps' in 
which water ice would be stable again s t sub lima- 
tion for billions of yea rs. IStacev et al. I (|l997l) and 
Campbell et al.l (| 1997T) found that the same signal 
would also result from a surface that was rough 
on the scale of the 13cm wavelength of the radi- 
ation and that even in sunlit regions this same 
high circular polarisation ratio was seen. Thus 
it is unclear whether or not water ice deposits 
are the explanation for these observations. Theo- 
ret ical mode l ling o f the radar scattering process 
by iFa et al. ( 2011 ) suggests that this technique 
will struggle to detect ice mixed in with lunar re- 
golith because of their simila r relat ive permittiv- 
ities. However, ISpudis et al. I (|2010h argue, based 
upon measurements from Chandrayaan-l's Mini- 
SAR experiment, that surface roughness cannot 
explain the signal seen from within some 'anoma- 
lous' polar craters where no signal is seen from 
the surrounding region, which might be expected 
to have had post-cratering roughness weathered 
away at a similar rate. 

An important method that has provided evi- 
dence relevant to the discussion of lunar water 
over the last 13 years is neutron spectroscopy. 
It was realised 50 years ago that the neutron 
flux leaking from the Moon encoded information 
about the hydrogen con tent in the lunar surface 
([Lingenfelter et ah 1961 ). Free neutrons are cre- 
ated via interactions between incoming cosmic 
rays and nuclei within the top few metres of the 
lunar regolith. Having been knocked out of nuclei, 
these neutrons subsequently evolve by scattering 
off other nuclei until they either escape through 
the surface or have sufficiently little kinetic en- 
ergy that they are absorbed by another nucleus. 
The fast neutrons (energies exceeding 0.5 MeV) 
predominantly lose energy through inelastic col- 
lisions that result in gamma rays being emitted 
from the temporarily excited nuclei that they hit. 
At intermediate energies (0.3 eV< E < 0.5 MeV), 
epithermal neutrons donate kinetic energy to their 
targets, which recoil following elastic collisions. As 



hydrogen nuclei have a similar mass to that of the 
neutron, they are particularly effective at draining 
energy rapidly from epithermal neutrons and send- 
ing them down to thermal energies (E < 0.3 eV), 
where they become increasingly susceptible to nu- 
clear capture. Thus, a hydrogen rich top metre of 
regolith would be expected to produce a deficit in 
the epithermal neutron flux leaking from the lunar 
surface. 

This first calculation of the anticipated sensitiv- 
ity of orbital neutron measurem ents to the regolith 
hydro gen c ontent was refined bv | Metzeer fc Drakel 
(|l990h and iFeldman etaH 4l99l|) with the bene- 
fit of improved knowledge of the regolith com- 
position. These models were confronted with 
observational evidence when results from Lunar 
Prospec tor's Neutron Spectr ometer were pro- 



duced (jFeldman et al.l Il998bl ). A 



convincing 

deficit of epithermal neutrons was detected in 
the vicinity of both lunar poles. However, Lunar 
Prospector (LP) was an omni-directional detector 
that accepted neutrons from all directions. Con- 
sequently, the spatial resolution of the published 
hydrogen maps was sufficiently poor that the ques- 
tion of how well these count rate dips correlated 
with surface features, such as craters, remained 
tantalisingly unanswered. Other results from 
LP's Neutron and Gamma Ray Spec t romete rs 
(lLawrence et all Il998t IFeldman et all Il998al |bl 
Elph ic et al.l 1998| ) provided a detailed and co- 
herent picture of the lunar surfac e composition. 



The analysis of the full data set (jFeldman et al 



2000, 120011 ). along with a detailed account o f the 
data reduction process ( Maurice et al. 2004 ) and 
improved modelling of both the d etector response 



and variation with composition (jLawrence et al 



2006) provide a benchmark against which current 



and future experiments can be compared. 

A recent development in the story of lunar wa- 
ter has come courtesy of measurements of the 
infra-red reflectance of sunlit regions of the Moon. 
Three different studies of the infra - red spectrum 
(IClarkl l2009t ISunshine et~aH 120091: IPieters et af 



2009) found evidence for absorption by oxygen 



to hydrogen chemical bonds on the lunar surface. 
Furthermore, the depth of the absorption features 
appeared to vary with solar zenith angle. Such a 
result would imply the existence of a contribution 
by hydroxyl and water molecules to the dynamic, 
if tenuous, lunar exosphere. This idea of a lunar 
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wat er cycle had been c o nside red by iButler (1997) 
and lCrider fc Vondrakl (|2000h . 

The most compelling evidence for the existence 
of lunar water came from the Lunar Crater Ob- 
servation and Sensing Satellite (LCROSS), which 
trailed a spent Centaur rocket into a permanently 
shaded part of the Cabeus crater near to the lu- 
nar south pole. Nine instruments on LCROSS 
analysed the plume of material resulting from this 
impact. It was determined to have contained 
water molecules at 5 .6 ± 2.9 per cent by mass 
(|Colaprete et al.ll2010f ). This uncertainty is domi- 
nated by the uncertainty in the total plume mass, 
and the detection of water vapour and water ice 
within the instrumental field of view was a much 
more significant 155 ± 12kg. 

Given that water does exist on the Moon, the 
questi on of how it ca me to be there becomes rel- 
evant. Arnold! ( 19791 ) considered various potential 
sources: delivery by comets or meteorites, solar 
wind reduction of iron in the lunar regolith or de- 
gassing from the lunar interior. As a result of the 
slow weathering of the lunar surface, the quantity, 
isotopic ratios and distribution of water, both with 
depth and position on the surface, should still con- 
tain information about the history of water deliv- 
ery to the Moon. It is thus important to make 
an accurate map of the distribution of lunar water 
in order to help disentangle the importance of the 
various delivery options. 

With this aim in mind, there are two differ- 
ent approaches that have been taken to try to im- 
prove upon the original LP hydrogen distribution 
maps. The software route involves using image 
reconstruction techniques to undo the instrumen- 
tal blurring by the omni-directional LP detector, 
whereas the hardware route is to make the detec- 
tor itself collimated. 



Elphic et al.l (|2007t) introduced a pixon-based 



image reconstruction algorithm that extracted 
more informa tion from the LP epithermal neutron 
data set, and lEke et al. I (120091) showed that, when 
all permanently shaded regions at both poles were 
stacked together, these data statistically favoured 
a reconstruction where the hydrogen was con- 
centrated into the permanently shaded regions. 
This result has important implications for the 
molecular form of the hydrogen, its method of de- 
livery to the Moon and survival on the surface. 
With improved maps of the permanently shaded 



polar regions from the Japanese SELENE mis- 
sion (|Noda et al.ll2008l ). lTeodoro et al.l (|2010h con- 
firmed this earlier result and made lunar polar hy- 
drogen maps that were given to the LCROSS team 
to help them decide where to cra sh their Centaur 



rocket in the search for w ater ice (|Colaprete et al 
20101 : ISchultz et al.ll2010h . 



The alternative approach, to deploy hardware 
that only accepts neutrons from a small fraction 
of 4-7T steradians, was taken with the Lunar Ex- 
ploration Neutron Dete ctor (LEND), described by 
iMitrofanov et all ( 20081 ) , on the Lunar Reconnais- 
sance Orbiter (LRO). Blocking epithermal neu- 
trons from other directions is a challenge best met 
by placing many moderating and absorbing nuclei 
in the way. However, the necessary requirement 
not to send large masses away from the Earth op- 
poses this straightforward approach to producing 
a collimated neutron detector. An attempt to bal- 
ance these competing constraints has been made 
with the 4 LEND Collimated Sensors for Epither- 
mal Neutrons (CSETN1-CSETN4), hereafter re- 
ferred to collectiv ely as CSETN. T his is one of 
seven instruments ( Chin et al"1l2007l ) on board the 
LRO mission that has, to date, been mapping the 
Moon since mid-2009. 

There has been some disagreement in the lit- 
erature concerning the s ensiti vity of the LEND 
CSETN. ILawrence et aD (|2010h argue that the de- 
tected flux should contain about 0.18 collimated 
neutrons per second. This conclusion was ar- 
rived at following both analytical calculations, us- 
ing geometrical considerations to scale the LP 
count rate, and Monte Carlo neutron transport 
modelling, which also provided an estimate of 
the uncollimated lunar background flux of 0.09 
counts per second. These findings were based 
upon the not entirely specific des cription of the de- 



tector configuration provided by IMitrofanov et al 
(|2008l ). The estimate for the collimated neutron 



count rate given by IMitrofanov et al. ( 2008 ) is 0.9 
counts per second. They also state that the total 
background from all sources should provide about 
1.2 neutrons per second of background counts. 
More recently, using data gathered by the LEND 
CSET N during the first year in orbit around the 
Moon, IMitrofanov et al" (l2010bl) assert that the 
collimated neutron count rate is more like 1.9 neu- 
trons per second, with 0.3 counts per second of un- 
collimated lunar signal and 2.8 neutrons per sec- 
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ond as a r esult of cosmic r a y strik es on LRO itself. 
However, lLawrence et al.1 (|201ll ) use the compo- 
sitional variation of the high- and low-energy ep- 
ithermal neutron count rates to estimate that the 
uncollimated lunar flux exceeds 2 counts per sec- 
ond, leaving very little collimated signal. 

The abundance and distribution of near-surface 
lunar hydrogen provides a valuable constraint on 
models describing the delivery and survival of 
volatile materials over billions of years. Making 
accurate maps of the lunar hydrogen distribution 
is thus one of the foremost aims for those seeking 
to understand this important aspect of the history 
of the inner solar system. As the ability to detect 
small dips in epithermal neutron flux is paramount 
when using orbital neutron data to infer hydrogen 
abundances in the top ~ lm of lunar regolith, it 
is crucial to know precisely what is the signal-to- 
background ratio for the LEND CSETN. The most 
pertinent question to address is what fraction of 
the detected neutrons actually originate from the 
lunar surface in the small field of view of the col- 
limator? The remainder of the detected neutrons 
would be uncollimated and originating either in 
the Moon or in the spacecraft itself. 

The Geosciences Node of NASA's Planetary 
Data System (PDS^) is a public repository that 
archives and distributes digital data related to the 
study of the Moon, and it contains data from 
the LEND CSETN. It is thus timely to compare 
these new data with those collected over a decade 
ago by the LP Neutron Spectrometer. The count 
rate variation with spacecraft position and alti- 
tude provides a direct means to decompose the 
total detected flux into collimated signal and un- 
collimated background noise. The main purpose 
of this paper is to answer the question, what frac- 
tion of the LEND CSETN counts are epithermal 
neutrons that truly originate from the small piece 
of lunar surface geometrically accessible to the col- 
limator? 

Section [2] of this paper provides a description 
of the data that will be used in this study and 
the likelihood analysis that is employed to deter- 
mine what fraction of the LEND CSETN count 
rate is collimated epithermal neutrons. Section [3] 
contains the results of using the data to decompose 
the total counts into those from the various com- 



ponents, as well as a discussion of how robust these 
findings are to potential shortcomings in the mod- 
elling of the count rate. These results are com- 
pared with the componen t fractions advocated by 
Mitrofanov et al.l ( 2010bl ). The findings are sum- 
marised in Section |4] 

2. Method 

This section contains a detailed description of 
the data being used, the model that is employed 
to describe them, and the procedure used to find 
the best-fitting model parameters. These param- 
eters are the fraction of the total LEND CSETN 
counts coming from each of the different possible 
sources, namely collimated lunar neutrons, uncol- 
limated lunar neutrons that have enough energy 
to pass through the sides of the collimator and 
uncollimated spacecraft neutrons. 

2.1. Data 

The data analysed in this paper are the first 
15 months of 'rectified science data' available 
within the LEND Reduced Data Record (RDR) 
data product available at the Geosciences Node of 
NASA's PDS. Data from all 4 LEND CSETN de- 
tectors are studied. The neutron count rate com- 
prises the sum of counts in channels 10-16, because 
these channels mostly contain neutron detections 
rather than the gamma rays or charged particles 
that prevail in channels 1-9 (Mitrofanov, private 
communication). To avoid needing to know de- 
tailed information concerning neutron transport 
within LRO and LEND, only the observations for 
which the nadir angle is less than 5° are used. 
In this configuration, the LEND instrument lies 
slightly beneath and to one side of the LRO space- 
craft and has an unobstructed view of the lunar 
surface. Other than on the first day of collection 
of these data (15th September, 2009), LRO was in 
a reasonably circular mapping orbit with altitudes 
generally between 40 km and 60 km. During the 
initial day of data, the spacecraft had yet to be 
put into its mapping orbit, so reached altitudes as 
low as 30 km and as high as 210 km. Throughout 
the rest of this paper, the zero of time will be 
defined as the start of the data collection period. 



1 http: / / pds-geosciences.wustl.edu 
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□ays since 15th Sept, 2009 



Fig. 1. — The box-car smoothed (with a box containing 10000 consecutive 1 second observations) LEND 
CSETN data for the first 15 months of operation. 
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Fig. 2. — A Mollweide projection of the raw 
LEND CSETN channel 10-16 data. This uses the 



HEALPix pixellation (jGorski et al.ll2005j ) at a res- 
olution of A^idc = 64, which corresponds to a pixel 
side length of ~ 28 km. 



2.1.1. Long-term temporal variations in the time 
series data 

The box-car smoothed time series data are 
shown in Fig. [TJ It is apparent that there are some 
glitches in these smoothed count rates. These 
are almost always associated with periods when 
no data are available from the detectors, and 
have a periodicity of roughly 2 weeks, or approx- 
imately half a lunar month. After the glitches, 
the count rate can be as much as 20% lower 
than before, and it gradually rises over a period 
as long as a few weeks. FigJ^] shows a map of 
the raw count rate. It is apparent that there 
are bands of low counts at —110° <longitude< 
-90° and 70° <longitude< 90°, and bands of 
high count rate at —80° <longitudc< —0° and 
100° <longitude< 180°. The detectors are orbit- 
ing in a plane whose normal is approximately the 
Earth-Moon axis when they are switched off. 

These artefacts are largely removed when the 
temporal variability is modelled as described in 
the next section. Recorded count rates are some- 
times unusually high both slightly before and after 
the glitch times, so data are removed in the peri- 
ods 5 minutes prior to the time when the data stop 
and 1 minute after the time when they reappear. 

The one example of a glitch where the count 
rate jumps upwards is found at t ~ 122.5 days, 
within a short period between two other glitches. 
To avoid any problems this might introduce, all 
data with [120 < i/days < 123] are discarded. 

Other notable features of Fig. Q] are that there 
are gaps without data at t ~ 80 and 250 days, 
after which the flux is significantly lower. As de- 
scribed on the PDS, this is because CSETN3 and 
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CSETN2 respectively are switched off at these 
times. Both of these detectors recommence record- 
ing data around t ~ 335 days. Along with the tem- 
poral variation described above, this should also 
be accounted for when trying to predict the total 
counts for any particular observation. 

The number of Is observations remaining after 
pruning the data as described above is ~ 35 mil- 
lion. They span a total of 456 days from 15th 
September, 2009. Various subsets of these data 
are used in Section [3] to investigate the robustness 
of the results, but the default in this paper is to 
use all of these observations. 

2.2. Model 

There are 3 components that will contribute to 
the LEND CSETN channel 10-16 count rate. The 
magnitudes of each component will depend upon a 
combination of the spacecraft latitude, longitude 
and altitude, a. The components are: 

1) collimated lunar epithermal neutrons, 

2) uncollimated lunar neutrons, with energies 
high enough to make it through the walls of the 
collimator, 

3) spacecraft-derived neutrons also entering into 
the detectors. 

The final component may also include any con- 
tribution to the channel 10-16 count rate due to 
charged particles resulting from cosmic rays strik- 
ing the spacecraft. 

The 3 He tubes used to count incoming neutrons 
in LEND, and also LP, are covered in cadmium to 
keep out low energy, thermal neutrons at energies 
below ~ 0.4eV. The e nergy dependence of both 
the detector efficiency ( Lawrence et al.l 20021) and 
the l unar neutron leakage flux ( McKinnev et al.l 
2006) means that low energy epithermal (LEE) 
neutrons with ~ 0.4 — lOOeV provide the bulk of 
the counts in these detectors. In addition to the 
Cd coating, the LEND CSETN are placed inside 
a collimator that consists of layers of polyethylene 
to moderate the neut rons and 10 B to absor b the 
moderated neutrons. Lawrence et al. esti- 
mate that the collimator is largely transparent for 
neutrons that strike it with at least ~ lOkeV. In- 
cident neutrons at these higher energies could be 
moderated to a low enough energy that the 3 He 
counter is able to detect them. It is this route 



that neutrons comprising components (2) and (3) 
must take in order to contribute to the LEND 
CSETN count rate. Component (2) reflects the 
high energy epithermal (HEE) neutron flux from 
the lunar surface; an interesting science product 
in itself. 

The technique employed here essentially con- 
structs a model for the total count rate that the 
LEND CSETN should measure in each Is obser- 
vation from each of the components. This model 
is based upon the LP maps of fast, moderated, 
epithermal and thermal neutron flux and a spa- 
tially invariant spacecraft background. Each of 
the model component fluxes also requires the in- 
clusion of both an altitude dependence and the 
temporal count rate gain variation, g(t). Once the 
individual components have their latitude, longi- 
tude, altitude and time variations computed, they 
are each normalised such that a fraction f\ of the 
total counts summed over all observations in the 
data set comes from component i. In this nomen- 
clature, determining /; is the purpose of this pa- 
per. The LEND CSETN data are used to deter- 
mine the relative probabilities of the different sets 
of component fractions, /i, i.e. they determine 
which component fractions provide the best match 
to the time series measurements. 

The evaluation of the spatial, altitude and tem- 
poral variations of the model component fluxes are 
described in the following subsections. 

2.2.1. Constructing the spatial variation in the 
model 

Given the well documented and well understood 
LP maps for the spatial variation of the thermal, 
epithermal and fast neutron fluxes, these data sets 
are the natural resource for designing whole- Moon 
templates of the various components of the LEND 
CSETN flux. Each component of the model re- 
quires a template that describes how the rela- 
tive count rate varies with longitude and latitude. 
These templates will hereafter be denoted with Ui 
for component i, and their normalisation is dis- 
cussed at the end of subsection 12.21 To model the 
spatial variation of the components of lunar ori- 
gin, various different resolution maps of the LP 
neutron data are made using the HEALPix pixel- 
lation scheme ( Gorski et al.ll2005l ). This partition- 
ing of the sphere has 12 equal area base pixels and 
finer pixellations are achieved by splitting each of 
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Fig. 3. — A Mollweide projection of the LP ep- 
ithermal neutron map including the pixon recon- 
structions in the polar regions, which is used as the 
template for component (1) of the model. This 
uses the HEALPix pixellation at a resolution of 
^Vgide = 64, which corresponds to a pixel side 
length of ~ 28 km. 
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Fig. 4. — A Mollweide projection of the thorium- 
corrected ep\ tho map, used as part of component 
(2) in the model. This uses the HEALPix pixella- 
tion at a resolution of A s jd = 64. 



the pixels into 4 equal area pixels. The level of 
resolution can be denoted by Aside, the number of 
pixels along a side of a base pixel. In this nomen- 
clature, the total number of pixels describing the 
spherical surface is A pix = 12A s 2 idc and, for the 
Moon, the typical pixel side length for the kth 
level of refinement is — (1778/A s id e ) km. 

Given that LRO can be at altitudes from 30 — 
210 km, although it spends most of the time at 
its nominal mapping altitude of 50 km, it is nec- 
essary to choose the appropriate resolution tem- 
plate, when calculating the model counts for a 
particular observation. This is done quite crudely 
by linearly interpolating between the pixel values 
either side of the FWHM for the component un- 
der consideration. There are inevitable inaccura- 
cies in this treatment, both because LP did not 
have comparable spatial resolution to component 
(1), and it did not produce a HEE map containing 
the same spectrum of neutron energies as would 
be expected to produce component (2). However, 
this simple approach will still provide an adequate 
comparison between LP and LEND CSETN data 
sets. 

It is necessary to find an expression for the vari- 
ation of the FWHM with altitude, so that an ap- 
propriate pixellation can be chosen to give a suit- 
ably smoothed measure of the lunar flux. This 
can be achieved by considering the neutron flux 
emanating from a small area on the lunar surface. 
If the detector is orbiting at an altitude a, a dis- 
tance r from the neutron source under consider- 
ation, with a corresponding arc length s between 
the sub-spacecraft point and the source, then the 
relative flux received can be written as 



/(*) 



/(* = 0) 



(^) 2 (cos0) 1+ «, (1) 



where a represents the effective beaming of neu- 
trons relative to the constant surface brightness, 
Lambertian case, and 9 is the angle from the 
surface normal to the spacecraft. This 'beam- 
ing' arises because the epithermal neutron num- 
ber density in creases with depth in the top layer 
of the regolith ( McKinnev et al.ll2006 ). and is anal- 
ogous to solar limb darkeni ng. Using Monte Carlo 
neutron transport models, Lawrence et al. (2006) 
found that a ~ 0.5 provided a good description 
of the beaming of epithermal neutrons from the 
lunar surface. Simple trigonometry yields the fol- 
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lowing expression for the change in flux with sub- 
spacecraft distance along the surface, from which 
the altitude dependence of the FWHM can be cal- 
culated: 



[ycos(s/r M ) - l] 1 " 



/(s = 0) V r M/ [y 2 + 1 - 2ycos(s/r M )]^ 



(2) 



where tm : 
dius and y 



1737.4 km represents the lunar ra- 
(Vm + a) I ru ■ This expression as- 
sumes that the detector has no preferred direc- 
tion. While not strictly true, comparison of the 
in ferred LP point - sprea d function (PSF) with that 
of Maurice et al.l ( 2004 ) suggests that it is an ad- 
equate description in that case. 

Having described how the component templates 
will be used, the next step involves actually cre- 
ating them. The template for the collimated LEE 
component should look like a noiseless, high spa- 
tial resolution version of the LP epithermal neu- 
tron map. While the spatial resolution of the LP 
epithermal map is that of an omni-directional de- 
tector, in the polar regions, the pix on reconstruc- 



tions of the LP epithermal data (jTeodoro et al 



2010) can be used to gain higher spatial resolution 
for the LEE component (1) template, u\, as shown 
in Fig. [31 Using the pixon maps allows the model 
to be extended from a resolution of -/V s id e = 64 to 
N sidc = 128 and 256 at |lat| > 71°. This highest 
resolution roughly corresponds to the anticipated 
footprint for component (1) in the model. The re- 
sults do not depend greatly on this detail in the 
model, as demonstrated in Section [3] 

The spatial variation of the HEE component 
(2) is somewhat less easy to estimate, both be- 
cause the relevant incident neutron energies de- 
pend upon precise details of the collimator, which 
are not available, and there does not exist a 
high-fidelity map of lunar neutron emission inte- 
grated with the appropriate energy-dependent fil- 
ter. However, neutron data from the LP Gamma- 
ray Spectrometer (G RS) anticoincidence shield 
(| Genet av et al. 2003) were used to derive infor- 
mation from epithermal neutrons in the energy 
range 100 — 500 keV, which has significant overlap 
with the LEND HEE neutrons. These 'moder- 
ated' neutrons showed a st rong correlation wit h 
LP-measured fast neutrons (jMaurice et al.l 120001) . 
which themse lves are correlated with average 
atomic mass ( Gasnault et al.l [2001). The mea- 



sured correlation of moderated neutrons with fast 
neutrons suggests that the HEE neutrons may 
also, in part, be correlated with fast neutrons. 
The LP moderated neutron data are consistent 
with neutron transport models that show HEE 
neutrons are highly corre lated with fast neutro ns 
and average atomic mass ( Lawrence et al.lfeoil ) in 
regions with low concentrations of hydrogen. For 
regions with enhanced hy drogen concentrations , 



modelled HEE neutrons (Lawrence et al. 2011) 



and measured moderated neutrons ( Genetav et al 



2003) show a decrease in neutron count rates that 
is very similar to LEE neutrons. 

Based on this information, a variety of educated 
guesses can be made for a HEE template by con- 
sidering all of the LP thermal, epithermal, mod- 
erated and fast neutron maps. In this study, two 
different types of initial HEE template have been 
devised. For the first t emplate, a map similar to 
the epi* constructed bv lFeldman et al. ( 1998b ) is 
created by subtracting 6.4% of the LP thermal 
count rate from the LP epithermal map, which 
effectively removes the LEE compositional varia- 
tion due to neutron absorption at energies below 
100 eV. To account for a fast neutron contribution 
as described above, a to-be-determined amount of 
the LP fast neutron map is added. The precise 
fraction of fast neutrons depends on details of the 
LEND construction, and is therefore left for the 
data to select. 

An option of correcting for residual, nearside 
neutron absorption in the epi* map is also con- 
sidered by subtractin g a portion of the LP tho- 
rium concentrations ( Lawrence et al. 2003h that 
spatially correlate with the residual neutron ab- 
sorptions. Fig.@]shows the template for this epi t/l ° 
component of the HEE variation and Fig. [5] shows 
the LP fast neutron map, which also contributes 
to U2, the template for component (2). 

The second method involves replacing all re- 
gions of the LP epithermal map where the LP 
thermal counts are lower than 20 neutrons per sec- 
ond with a scaled version of the LP fast neutron 
map. This is essentially targetting the mare re- 
gions, where the LEE and thermal neutrons have 
a dip in count rate, whereas the fast neutrons have 
an excess as would be expected from LP mod- 
erated neutron measurements and HEE neutron 
transport simulations. Unlike for component (1), 
the HEE lunar template only needs to be con- 
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structed at the low resolution corresponding to an 
uncollimated detector. 

By using only nadir facing observations, the 
spacecraft background should be independent of 
location over the surface, although not indepen- 
dent of altitude, so u 3 is taken to be a constant, 
independent of latitude and longitude. 




2.2.2. Constructing the altitude dependence in 
the model 

The altitude dependence of the lunar flux com- 
ponents, h\ and hi, can be found by integrating 
eq. ([2]) over the appropriate fraction of the lunar 
surface out to the horizon for the chosen altitude, 
a. For the collimated epithermal component, the 
limit of the integration surface is approximated as 
s = a/3 co \, where /3 co i = 5.6° is the half-opening an- 
gle of the collimator. There is almost no altitude 
variation of the collimated component because the 
flux from each visible patch of surface drops like 
a -2 , while the area of surface grows like a 2 . As the 
beaming of neutrons is weak, the Moon presents a 
disc of almost uniform neutron surface brightness, 
these two factors almost entirely cancel out, pro- 
vided that the Moon always fills the collimator's 
field of view. 

The uncollimated lunar HEE component is 
computed in a similar way by integrating out to 
the horizon. In this the spacecraft altitude 

increases, the Moon subtends a smaller angle and 
the uncollimated flux decreases. One might expect 
that the beaming of HEE would differ from that 
of LEE, because HEE neutrons are more likely 
to have been created by cosmic rays hitting the 
lunar surface at more grazing angles. The model 
is given the freedom to choose a value of a ^ 0.5 
for component (2), if the data prefer. 

The altitude dependence of the spacecraft flux 
arises as a result of the solid angle of cosmic ray- 
containing sky blocked out by the Moon. Conse- 
quently, higher altitudes lead to more cosmic rays 
hitting the spacecraft and more locally generated 
neutrons being detected by the LEND CSETN. 
The flux of cosmic rays hitting the spacecraft as a 
whole varies like 



CRT 



1 



1 
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(3) 
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13 14 15 16 

LP fast neutron count rate/s" 



Fig. 5. — A Mollweide projection of the LP fast 
neutron map, used as part of component (2) in 
the model. This uses the HEALPix pixellation at 
a resolution of A s id c = 64. 




One might wonder if the flux of cosmic rays hitting 
the nadir facing side of LRO is more relevant for 



100 150 
AHitude/km 

Fig. 6. — The altitude dependence, h(a), of the 
count rates for the various different components. 
A solid line traces the variation with altitude of the 
collimated LEE lunar flux. Uncollimated HEE lu- 
nar flux, h 2 , drops off as shown by the dotted line. 
The spacecraft flux, h 3 , is traced by the dashed 
line. 
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predicting the neutron flux reaching LEND. This 
varies with altitude more rapidly: 



CRF 



1 

y 



(4) 



An experiment was performed where two compo- 
nents of uniform background were introduced, one 
with no altitude dependence and another varying 
like eq. (|4]), these being the extreme limits of how 
component (3) might depend upon altitude. The 
data picked out a most likely model with a combi- 
nation of these two components that corresponded 
to an overall spacecraft background altitude de- 
pendence as given by eq. ([3]). Thus, the altitude 
dependence of component (3) was set to 



h 3 oc 1 



1-- 

y 



(5) 



Figj6] shows the altitude dependence of all three 
model components of the count rate, with each 
normalised to unity at a = 50 km. 

2.2.3. Constructing the temporal variation in the 
model 

The variation of the effective count rate gain of 
the detector must be accounted for in the model 
to describe the received flux as a function of lon- 
gitude, latitude, altitude and time. While the rea- 
son for the glitches and gradual recoveries in the 
count rate is unclear, physical systems might be 
expected to have exponentially time-varying be- 
haviour. Thus, a natural and effective function to 
fit to each interglitch period is 



C = C + A{1 - e 



_ „-A(t-to) 



(6) 



where C represents the count rate, t, the time, Co 
is the count rate at a reference time to m this in- 
terglitch period, and A and A are constants to be 
determined by fitting to the data. This approxi- 
mate correction is applied for the first evaluation 
of the most likely component fractions. It takes 
into account the exponential recoveries in count 
rate and the sharp drops across periods when the 
detectors are not returning measurements. If there 
are any cosmic ray flux variations on timescales of 
a week or longer, then they will also be absorbed 
into this empirical temporal correction. 

Any temporally correlated residual differences, 
on timescales shorter than about a week, between 



the most likely model and the data are assumed 
to result from cosmic ray flux variations. These 
are estimated using the boxcar smoothed data and 
model over periods of 40ks, and this 'cosmic ray' 
factor can also be placed into the temporal vari- 
ation in the model. It is, of course, possible that 
inexact glitch removal can lead to incorrectly in- 
ferred cosmic ray variation, as these two factors 
are essentially degenerate. Only the assumption 
that the glitches are described by e q. (JSJ) decouples 
the two. iMitrofanov et all (|2010bl) do not mention 
the glitches in the data despite the fact that they 
require systematic corrections more than an order 
of magnitude greater than the quoted statistical 
uncertainties in their results. A comparison of the 
inferred relative cosmic ray flux with an indepen- 
dent proxy, namely the terrestrial neutron flux, is 
given in Section [3] 

To take into account the combination of the 
long-term variations in the data, count rates for 
each model component are normalised to produce 
the observed total number of counts over the en- 
tire data set when /; = 1. This corresponds to 
3.58 counts/second. As a result, each observation 
time has its own normalisation factor, defined such 
that the measured count rate is a factor g(t) times 
what it would have been without these temporal 
variations associated with the detectors. For the 
initial likelihood evaluation, g(t) only corrects for 
the varying number of detectors and the glitches in 
the data, whereas for the second likelihood eval- 
uation it also includes the cosmic ray correction 
factor. 

Taking into account all of these different com- 
ponents of the modelled count rate, the antici- 
pated count rate for any given observation, j, can 
be written as 



(7) 



The free parameters are the component fractions, 
/i, determining the relative contributions of the 
components. In fact only 2 of these are indepen- 
dent because the sum of all fractions is 1 in order 
to recover the correct mean neutron count rate. 

2.3. Likelihood evaluation 

By choosing a particular set of fractions, {fx}, it 
is possible to determine what is the probability of 
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measuring the time series LEND CSETN data set. 
Varying {/;} allows the most likely set of compo- 
nent fractions to be inferred. The probability of 
measuring dj counts in observation j, where the 
model predicts mj counts, is given by the Poisson 
distribution: 



If the likelihood of measuring a set of n bs data 
values is denoted by 

"obs 

L =H p v (9) 

j=i 

then, ignoring the model-independent denomina- 
tor, a log-likelihood can be defined as 

"obs 

lnL = ^(d j ln(m j )-m j ). (10) 
j=i 

The log-likelihood is calculated throughout the 
relevant region of the component fraction hyper- 
space, and the most likely parameter values, as 
defined by the data, are located. These log- 
likelihoods can be converted back to probabilities 
that the LEND CSETN data would arise given this 
model, and normalised so that the probability in- 
tegrated over parameter space equals 1. The one- 
dimensional marginalised probability can be cal- 
culated by integrating these probabilities at fixed 
values of any particular parameter, An addi- 
tional constraint is that none of the components is 
allowed to produce a negative contribution to the 
count rate. 

3. Results 

This section contains the results of applying 
the method described in Section [5] to the LEND 
CSETN data. Before describing the inferred com- 
ponent fractions, two additional improvements 
to the previously outlined model are described. 
These are essentially motivated by information in 
the data themselves. 

3.1. Additional empirical corrections 

Using the HEE template based upon the com- 
bination of the thermal-corrected LP epithermal 
and LP fast neutron maps, the entire 15 months 



of glitch-corrected LEND CSETN were used to de- 
termine the maximum likelihood component frac- 
tions. The temporal correlations in the residuals 
between the data and most likely model, on time 
scales of 40ks, are assumed to arise as a result 
of variations in the cosmic ray flux. This implies 
that, with the residual for each time series mea- 
surement defined as 

Rj=dj-mj, (11) 

periods when the residuals are systematically pos- 
itive or negative correspond to times when the cos- 
mic ray flux is high or low respectively. 

A test of this assumption is shown in Fig.[7]with 
the correlation of the ratio of the smoothed LEND 
CSETN data to the smoothed model values and 
the neutron flux measured on Earth at the Thulc 
station operated by the Bartol Research Institute. 
It is apparent in the time series in Fig. [7] that high 
neutron fluxes are often coincident with periods 
where the LEND CSETN data are systematically 
higher than the most likely model predicts. Con- 
sidering the ratio of smoothed data to smoothed 
model as a proxy for cosmic ray flux, the rela- 
tive mean values for each day can be compared 
with the ground-based data. A straightforward 
test gives a Spearman rank correlation coefficient 
of 0.19, which is highly significant. 

Fig. [7] also shows a long-term trend in the Thule 
data that is not present in the LEND CSETN 
cosmic ray proxy. This is because that temporal 
variation has been absorbed into the glitch correc- 
tion described earlier. Consequently, a fairer com- 
parison results from normalising the ground-based 
count rates by their one-month box-car smoothed 
values. This effectively removes the long timescale 
trends from the Thule data. In this case, the 
comparison between the LEND CSETN cosmic 
ray proxy and the detrended Thule neutron count 
rate gives a Spearman coefficient of 0.38. A com- 
parable test for data from other monitoring sta- 
tions gives Spearman coefficients of 0.32 (Swarth- 
more/Newark), 0.29 (McMurdo) and 0.27 (South 
Pole), all of which are highly statistically signifi- 
cant. 

This evidence provides an important consis- 
tency check for the model used here - namely 
that one can infer something like the cosmic ray 
flux from just the LEND CSETN count rates and 
the approximate glitch corrections represented by 
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1 day— averaged daia/mode 
i day- averaged Thule neutron count raLe 



200 s50 300 
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Fig. 7. — Time series showing the variability in the 1-day averaged neutron flux at the Thule station of the 
Bartol Research Institute and the ratio of the 40ks box-car smoothed LEND CSETN data to the similarly 
smoothed model. Both of these quantities can be viewed as proxies for the variability in the cosmic ray flux. 



eq. ([6]). Using the temporal correlations in the ra- 
tio of smoothed data to smoothed model on scales 
of 40ks as a proxy for the variability of the cos- 
mic ray flux, the correction factor, <?(ij), is mul- 
tiplied by this cosmic ray correction factor to re- 
move these correlated residuals, and the maximum 
likelihood component fractions are recalculated. 

This second iteration makes very little differ- 
ence to the inferred component fractions, although 
it does significantly improve the match of the most 
likely model to the time series data. A map of the 
residuals from this second iteration is shown in 
Fig. [HI The pixel size is ~ 220 km, corresponding 
to AT side = 8 in the HEALPix notation. While the 
temporal correlations in the residuals have been 
removed, there are still significant spatially corre- 
lated residuals. The collimated flux should have a 
spatial dependence that is well described by com- 
ponent (1) in the model, but, as noted in Section[2] 
the HEE component is more of a challenge. The 
most likely reason for these residuals is an inap- 
propriate choice of model HEE template, so it is 
updated to remove all of the residuals measured 
on the large scales shown in Fig. [8] The largest 
correction factor in the template count rate aris- 
ing from this change is only ~ 3%. While this de- 
creases the number of degrees of freedom provided 
by the data, it does not preclude the residuals on 
smaller scales subsequently being used to test how 
well the model describes the data. 

In contrast to the temporal correction, which 
could be cross-checked against an independent es- 
timate of the cosmic ray flux, there is no other 
available map of the Moon with the appropri- 




-4-3 3 4 

r=(d-m)/ a 

Fig. 8. — The spatial correlation of the reduced 
residuals (residuals normalised by an uncertainty 
in the pixel, tr, derived from Poisson counting 
statistics) for the default HEE template, having 
only corrected for the cosmic ray flux variations. 
Blue regions represent areas where the model over- 
estimates the LEND CSETN data. N sidc = 8 is 
used for this map, corresponding to pixel scales of 
- 220 km. 
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ate energy-dependent filter applied, whatever that 
might be in detail. One reason to think that this is 
not greatly biasing the results of this study is that, 
when starting from the alternative HEE template 
described in the previous section, the residuals are 
significantly improved relative to those in Fig. [8] 
and despite this, the inferred component fractions 
hardly change. Thus the results found here are 
robust to choosing very different input HEE tem- 
plates. 



200 



150 



. 100 



50 



j a 



0.2 



0.4 



0.6 



Fig. 9. — Marginalised probability distributions 
for the fractions of the total LEND CSETN counts 
in components: 1 - collimated lunar LEE neu- 
trons (red), 2 - uncollimated lunar HEE neu- 
trons (blue), and 3 - locally produced spacecraft 
neutrons (black). The coloured vertical dotted 
lines are the co r respon ding values asserted by 
Mitrofanov et all (|2010bl ). 



3.2. Most likely component fractions 

Fig. IH1 shows the marginalised ID probability 
distributions for the component fractions after 
both the cosmic ray correction and the HEE map 
modification. The most likely component fractions 
only change from {/i,/ 2 ,/ 3 } = {0.02,0.44,0.54} 
to {0.01, 0.445, 0.535} as a result of both of these 
changes. Statistical uncertainties on these values, 
as represented by the widths of the probability dis- 
tributions are barely more than ~ 0.01. There are 
larger systematic uncertainties that are considered 
in more detail later in this secti on. The preferred 
compon ent fractions reported bv lMitrofanov et al 



(|2010bl ) are also shown in this figure. Including 
the range of values they quote, their / 3 is consis- 
tent with that found here, namely just over half of 
the LEND CSETN count rate comes from neu- 
trons originating in the spacecraft itself. How- 
ever, the division of the lunar flux between com- 
ponents (1) and (2) is e ffectively reversed, with 
Mitro fanov et al. (|2010bh claiming that the colli- 
mated flux forms the vast majority, whereas the 
LEND CSETN data, coupled with the simple 
model based on previous knowledge from LP pre- 
sented in this paper, suggests that the lunar flux 
has roughly ten times as much coming through the 
collimator walls in a dominant component (2) as 
it does in collimated neutrons. 

One might wonder if the reason for this dis- 
crepancy is that the LEND CSETN detector it- 
self perhaps has a very different behaviour to the 
LP cpithermal neutron detector, despite the fact 
that they are both 3 He-filled Cd-covered tubes. If 
this were the case, then the model for component 
(1) could be in error, leading to incorrect com- 
ponent fractions. However, if one considers the 
map of lunar neutrons produced by the LEND 
SETN, which is an equivalent detector strapped 
to the outside of the LEND collimator, this expla- 
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Fig. 12. — Histograms of the reduced residuals, r, 
in pixels at a resolution of iVside = 32. The black 
curve is a Gaussian with width 1. Blue, cyan and 
green curves show the results for the most likely 
model before any corrections, after the cosmic ray 
correction, and after both cosmic ray and HEE 
template corrections respectively. The reduced x 2 
values, accounting for the decrease in the number 
of degrees of freedom caused by the data-based 
corrections, are given in corresponding colours in 
the legend. Results in red and m agenta are for the 
compon ent f ractions favoured b y Mitrofanov et al.l 
(|2010ri ) and lMitrofanov et all (|201ll ) respectively. 



Fig. 11. — Mollweide projection of the LEND 
SETN corrected count rate map, using iV S idc = 32. 
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Fig. 13. — The relative count rate as a func- 
tion of altitude. Points show the corrected LEND 
CSETN data (uncertainties are smaller than the 
symbol size). The blue, red and magenta lines 
show the trends from the most likely model found 
here and the as s erted c ompo nent fractions of 
Mitrofanov et all (|2010bl ) and iMitrofanov et al 
(|201 lh respectively. 



nation can be rejected. Figures [10] and [IT] show 
the corrected (for glitches and cosmic ray varia- 
tions) maps for the LEND CSETN and SETN re- 
spectively. It is clear from this comparison that 
the CSETN map is not a sharp version of that 
from the SETN, which is what would be produced 
if the lunar variability in the CSETN map were 
dominated by collimated neutrons. The lunar de- 
pendence of the CSETN map is driven by some- 
thing different to the SETN results. This is not 
co nsistent with the a sserted component fractions 
of IMitrofanov et al.l (|2010bi ) . Instead it reflects 
the preponderance of higher e nergy neutrons in 



the CSETN lunar flux. In fact, Mitrofanov et al 



(|201ll ) state that a significant fraction of their 'lo- 
cal background from the LRO spacecraft' is actu- 
ally scattered lunar HEE, which will have both a 
different spatial and altitude dependence to back- 
ground arising from cosmic rays striking LRO. 
They give {/i,/ 2 ,/ 3 } = {0.32,0.22,0.46} as the 
component fractions. 

One can study the distribution of pixel reduced 
residuals, r, (i.e. residuals normalised by the Pois- 
son counting error in each pixel) on scales smaller 
than that used for updating the HEE model tem- 
plate to see how well the most likely models fit 
the LEND data. Fig. [12] shows histograms of the 
pixel residuals for iV s ide = 32 for the most likely 
models before either cosmic ray or HEE template 
corrections, after each of these model upd ates, 
and for both the Mitrofanov et al.l (|2010bl ) and 
Mitrofanov et al. (|2011 ) component fractions fed 



into this corrected model. The most likely model 
improves significantly as a result of the tempo- 
ral and spatial corrections. However, even before 
cither of them, it was still vast ly superior to ei 



ther th e c omponent fraction s of | Mitrofanov et al 



(|2010bD or IMitrofanov et all (|201ll) . for which the 
jump in the distribution at r — 4 contains all pix- 
els with residuals at least this large. These are 
essentially pixels in the mare regions where the 
data show a high count rate while the collimated 
flux should show a decrease relative to other equa- 
torial locations. The number of degrees of freedom 
is such that, for a good fit, the reduced \ 2 should 
be lztO.Ol. Thus, even after both corrections, the 
most likely model found here is still not quite a 
'good' fit. However, given the approximate treat- 
ment of the spatial component of the modelling 
on scales of the instrumental PSF, this is an im- 
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pressive agreement considering the small number 
of model parameters being fitted to the data. 

Given the significant difference between the 
most likely compo nent fractions found here , and 
those advocated b y Mitrofanov et al. ( 2010bl ) and 
iMitrofanov et all ( 2011 ). one might expect that 
there would be other tests that would allow the 
data to choose between these two models. In par- 
ticular, given that the uncollimated lunar flux de- 
creases strongly with increasing detector altitude, 
whereas the collimated flux is approximately con- 
stant, the altitude dependence of the data should 
also be a powerful discriminant. There are cor- 
relations of typical spacecraft altitude with po- 
sition over the surface, such that LRO is usu- 
ally lower over the mare regions than elsewhere, 
so one needs to account for these by compar- 
ing each time series observation with the model 
count rate for that observation. These results, 
normalised by the mean overall count rate and 
binned with respect to altitude, are shown in 
Fig. [13] The data show a slight decrease of count 
rate with increasing altitude that can only be re- 
produced if the model contains enough uncolli- 
mated lunar flux to overcome the opposite trend 
brought from the dominant spacecraft background 
(see Fig. [6]). Once a gain, the asserted compo- 
nen t fractions of bot h Mitrofanov et al. ( 2010bl ) 
and lMitrofanov et al.l(|201lh are completely incon- 
sistent with the data. They imply an altitude de- 
pendence that is the opposite of that seen in the 
data. This points to the lunar flux predominantly 
being uncollimated rather than collimated. 

3.3. Investigation of potential systematic 
errors 

Given that the statistical uncertainties are so 
small, it is pertinent to ask how sensitive are these 
results to different assumptions underpinning the 
model. Table 1 contains results for various differ- 
ent model assumptions and subsets of the data, 
from which it is possible to gain some insight into 
the robustness of the inferred component fractions. 

The effect of not correcting the HEE template 
using the residuals is shown in row 2, whereas ne- 
glecting the 'cosmic ray' temporal correction too 
leads to the fractions in row 3. While these cor- 
rections greatly improve the residuals, as shown 
in Fig. I12[ they have no significant impact on 
the most likely component fractions. Rows 4 and 



5 show the inferred most likely component frac- 
tions for the cases where the initial HEE template 
does not include the thorium correction and for 
the alternative LP epi/fast template respectively. 
Once again, the component fractions are robust 
to these changes to the model. The same story is 
repeated for rows 6 and 7, where the decoupled 
pixon reco nstructions of t he LP epithermal map 

Of) are rep l aced w ith 




the coupled ones produced bv lEke et al.l ( 20091 ). or 
just the straightforward binned LP data. The fi- 
nal aspect of the model that was varied was the 
minimum pixel size used to bin the LP data to 
make the templates for the HEE model. Rather 
than using the default iV S ido = 64, one tier higher 
HEALPix resolution was adopted to model the 
HEE component. Once again, this made little dif- 
ference, despite gaining some resolution at the cost 
of introducing more noise in the model template 
for the HEE component. Changing the pixel size 
for the LEE component also had little effect, be- 
cause it is such a sub-dominant component. 

These results all provide evidence to support 
that none of the more uncertain aspects of the 
model are greatly impacting upon the component 
fractions inferred here. Keeping the default model 
choices, the data themselves were then split into 
different subsets, to see if similar results were ob- 
tained. 

Splitting the time series data into the 4 dis- 
tinct periods within which the number of active 
collimators does not change yields the results in 
rows 9-12. There seems to be a slight preference 
for data taken after t = 90 days to have a higher 
fraction in the collimated component, mostly at 
the expense of the spacecraft background. It is 
difficult to ascertain what aspect of the experi- 
ment altered to give rise to this change, which is 
statistically significant. 

Restricting the data to include only observa- 
tions made from restricted ranges in altitude does 
not appear to make any strong systematic changes 
to the inferred component fractions (rows 13-15). 
Using only the data for either the near or far side of 
the Moon yields almost identical most likely com- 
ponent fractions (rows 16 and 17), and the results 
for polar and equatorial regions in rows 18 and 19 
are very similar too. 

A significant difference is seen between the com- 
ponent fractions inferred from northern and south- 
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Table 1: Most likely parameter values calculated on a grid with spacing 0.005. Statistical uncertainties are 
of order 0.01 in all cases, which is smaller than the systematic variations arising from different subsets of the 
data, or model assumptions. The default model is for the 'thorium-corrected' input map. Lines 1-8 represent 
changes to the model, whereas lines 9-21 are for the default model with different subsets of the data. 
Row number Experiment /i fa 



1 


default 


0.010 


0.455 





.535 


2 


no HEE correction 


0.020 


0.450 





.530 


3 


no CR correction too 


0.020 


0.440 


0. 


.540 


4 


no thorium 


0.005 


0.460 





.535 


5 


epi/fast HEE template 


0.010 


0.465 





.525 


6 


coupled reconstructions 


0.010 


0.455 





.535 


7 


no pixon reconstructions 


0.005 


0.455 





.540 


8 


max iVsidc = 128 


0.010 


0.455 





.535 


9 


< t/days< 90 


0.000 


0.450 





.550 


10 


90 < i/days< 250 


0.040 


0.435 





.525 


11 


250 < i/days< 335 


0.040 


0.450 





.510 


12 


335 < i/days< 456 


0.035 


0.455 


0. 


.510 


13 


a/km< 45 


0.005 


0.445 





.550 


14 


45 < a/km< 55 


0.020 


0.450 





.530 


15 


55 < a/km 


0.015 


0.445 





.540 


16 


near side 


0.010 


0.460 


0. 


.530 


17 


far side 


0.005 


0.455 





.540 


18 


| latitude | > 60° 


0.020 


0.430 





.550 


19 


latitude < 60° 


0.015 


0.440 





.545 


20 


latitude> 


0.040 


0.415 





.545 


21 


latitude< 


0.010 
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ern hemisphere data (rows 20-21), with the north 
preferring more collimated and less uncollimated 
lunar flux than the south. This presumably re- 
flects residual inaccuracies in the HEE template. 
It is reassuring though that, while statistically 
significant, these systematic changes in the most 
likely component fractions are of a magnitude that 
is very small relative to the diffe rence between the 



values inferred here and those of Mitrofanov et al 



poiObl) 



3.4. Comparison with previous studies 

Breaking down the average ~ 3.6 counts per 
second (over the whole period used) to count rates 
in the various components, the results presented 
here suggest that on average, these comprise ~ 0.1 
per second in collimated epithermals, ~ 1.6 in 
uncollimated lunar HEE, and ~ 1.9 per second 
in spacecraft generated counts. Taking into ac- 
count the variable number of sensors active during 
the period, the collimate d epithermal count rat e 



is similar to the result of Lawrence et al. ( 2010l ) 



This is not surprising, because it is based on a 
simple geometrical argument that scales the flux 
from the omni-directional LP experiment to the 
solid angle subtended by the collimator. 

I n the supporting online m aterial for their pa- 



per, iMitrofanov et al.l (|2010bf ) perform a similar 
calculation where they scale the count rate from 
the LEND SETN to estimate the collimated lunar 
flux reaching the CSETN. However, they neglect 
the contribution of the locally generated back- 
ground to the SETN count rate and overestimate 
the beaming of lunar epithermal neutrons. Both 
of these factors will lead to an erroneously high 
anticipated collimated flux. 

The neutron beaming can be assessed using the 
type of argument that led to eq. ((2]). In order 
to recove r the beaming factor of 'about 2' as- 
serted by IMitrofanov et all (|2010bl ). the neutrons 
would need to be beamed with a ~ 1.5, in con- 
trast to the a = 0.5 measured from Monte Carlo 
neutron transport mo delling of the lunar regolith 



( Lawrence et al. 2006). This factor of 'about 2' 



appears to be 50% larger than one would infer for 
a = 0.5. 

Without an accurate description of the LEND 
instrument and detailed modelling of neutron 
transport within the spacecraft, it is difficult to 



quantify the extent of the locally-generated back- 
ground. However, it should be noted that the 
SETN is placed next to a large piece of polyethy- 
lene, which is known t o moderate and sca tter neu- 
trons very effectively ( Elphic et al. 20081 ). This is 
the reason it is used as the outer part of the col- 
limator. Thus it is implausible that there are not 
significant numbers of neutrons entering the SETN 
that come from this polyethylene, rendering any 
attempt to scale the SETN count rate to predict 
the CSETN count rate of little worth. Compar- 
ing the LEND SETN and LP epithermal maps in 
Figures [11] and [3] respectively, it is apparent that 
they look qualitatively different in a way that sug- 
gests the SETN has an extra HEE component that 
decreases the dip in the mare regions. 

While the collimated flux found here is in agree- 



ment with that anticipated by lLawrence et al 
( 2010h . there are approximately 20 times as 
many background lunar HEE neutrons than 
they predicted based on an incomplete knowl- 
edge of the dete ctor geometry. This appears 
to be a result of Lawrence et al. ( 2010h having 
been too optimistic in their assumption as to 
how effective t he col limator was. Specifically, 
(|201dh assumed a 10 B density of 



Lawrence et al 



2.5g/cm 3 , whi le the true density was s lightly less 



than lg/cm 3 dMitro 
tion, lLawrence et al. 



anov et al. 2010a). In add! 



( 2010h assumed that except 
for the opening aperture, the collimating mate- 
rial completely enclosed the 3 He sensors. In con- 
trast, the true LEND collimator has "holes" that 
allow neutrons to enter the 3 H e sensor through 



the back side of the collimator (Mitrof anov et al 



2010aJ). Both of these effects will significantly in- 
crease the true background comp ared to that es- 
timated by lLawrence et al" ( 2010l ). Consequently, 
one should view their conclusions concerning the 
LEND CSETN signal-to-noise ratios as having 
been too favourable by a factor of a/20 ~ 4. 

4. Conclusions 

A comprehensive maximum likelihood analysis 
has been performed on the LEND CSETN data 
to determine the sources of the neutrons that it 
detects. The fraction of detected neutrons that 
are collimated lunar epithermal neutrons is under 
0.05, with the most likely fraction being ~ 0.01. 
0.4 — 0.45 are uncollimated lunar HEE neutrons, 
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and the remaining 0.5 — 0.55 are neutrons gener- 
ated in LRO by direct cosmic ray strikes. The two 
lunar component fracti ons are almost the oppo - 
site of those asserted bv lMitrofanov et all ( 2010bf ). 
who favour the lunar flux being predominantly col- 
limatcd. However, the data themselves are com- 
pletely inconsistent with this interpretation, both 
using the full likelihood analysis and also from a 
simple inspection of the global S ETN and CSETN 
maps. The conclusions drawn bv lMitrofanov et al 



(|2010bf ) concerning hydrogen distributions should 
thus be viewed with skepticism. 

Given that the LEND CSETN is effectively 
producing an uncollimated map of the lunar HEE 
neutrons, albeit with a small contamination from 
a LEE collimated component and a large, but spa- 
tially uniform spacecraft background, the exper- 
iment would more usefully be called the LEND 
Sensor for High Energy Epithermal Neutrons 
(SHEEN). One would expect that the HEE neu- 
trons will be less beamed than the LEE ones, be- 
cause they are typically created by more grazing- 
incidence cosmic rays. This, in combination with 
the 50 km orbit of LRO relative to 30 km for LP, 
counteracts the effect of the small fraction of col- 
limated LEE neutrons, meaning that the effective 
FWHM for neutrons entering the LEND SHEEN 
is comparable with that of the LP omni-directional 
epithermal detector. 

The best available map of the variation of the 
lunar HEE flux, determined in this study and 
shown in Fig. 1 141 differs from the input tem- 
plate for the lunar HEE component that was 
estimated from LP results. This implies that 
the LEND SHEEN data contain extra informa- 
tion that, in conjunction with Monte Carlo neu- 
tron transport simulations, can be used to con- 
struct a global map of the lunar hydrogen distri- 
bution (albeit at a broad spatial resolution like 
LP). An ambiguity in the inferred hydrogen abun- 
dance from LEE neutron measurements in non- 



polar regions (jMaurice et al.ll2004l : lLawrence et al 



2006 ) will be absent from the HEE neutron case 
( Johnson et al. 20021 ). Specifically, HEE neutrons 
are, unlike LEE neutrons, unaffected by absorp- 
tion ( Lawrence et al.l 120061 ), but are still sensi- 
tive to hydrogen var iations ( Genetav et al.l 2003 ; 
Lawrence et al.ll201ll ). 

A global map of the hydrogen abundance would 
provide a valuable, new constraint on models for 



the lunar water cycle at low latitudes and a use- 
ful consistency check of the high-altitude hydrogen 
maps from LP. With the recent infra-red measure- 
ments indicating the existence of surficial water at 
low latitudes and the LCROSS discovery of water 
in a small patch of the Cabeus crater near to the 
lunar south pole, a global hydrogen map would 
help to tie together these two separate strands of 
evidence concerning lunar water. This study has 
shown that data from the LEND experiment have 
the capacity to provide a global hydrogen map, 
at a resolution similar to that which LP achieved 
near to the poles. However, this conclusion relies 
upon a correct interpretation of the observational 
data and understanding of the instrumental ca- 
pabilities. In order to exploit the full potential 
of the LEND SHEEN data with confidence, one 
would need a more detailed modelling of the de- 
tector response, with both position and energy, 
than is currently available. 
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Relative HEE count rate/s 1 



Fig. 14. — Mollweide projection of the final rela- 
tive HEE count rate map from the LEND SHEEN 
data set using N s - 1( i c = 64. The arbitrary normal- 
isation is chosen to match the mean count rate 
from the HEE component over the first 456 days 
of lunar mapping. Note that the correction for 
the correlated residuals on scales of ~ 200 km is 
visible. 
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